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Abstract 

In this work, we show that a large class of models with a composite dark sector undergo 
a strong first order phase transition in the early universe, which could lead to a detectable 
gravitational wave signal. We summarise the basic conditions for a strong first order phase 
transition for SU(N) dark sectors with Uf flavours, calculate the gravitational wave spectrum 
and show that, depending on the dark confinement scale, it can be detected at eLISA or 
in pulsar timing array experiments. The gravitational wave signal provides a unique test 
of the gravitational interactions of a dark sector, and we discuss the complementarity with 
conventional searches for new dark sectors. The discussion includes Twin Higgs and SIMP 
models as well as symmetric and asymmetric composite dark matter scenarios. 


1 Introduction 

Violent phenomena in the early universe can lead to large anisotropic fluctuations in the energy 
momentum tensor, which can act as sources for gravitational waves (GW). Strong first order 
phase transitions (PT) are an example of such a phenomenon, and it is well known that they can 
produce GWs [1-4]. Once produced, GWs propagate through space almost undisturbed, and can 
therefore serve as a unique probe of phenomena in the early universe. 

Phase transitions in particle physics are usually associated with symmetry breaking, i.e. 
with the transition of the universe from a symmetric phase to a phase of broken symmetry, as 
the temperature drops below a critical temperature T*. Within the standard model (SM) of 
particle physics, at least two phase transitions should take place, associated with the breaking of 
electroweak symmetry around T* ~ 100 GeV and with the breaking of chiral symmetry at the 
time of the QGD phase transition, T* ~ 0.1 GeV. 

Today we know that both the QGD and the electroweak phase transition are not first order, 
but proceed through a smooth cross-over [5-8] , and can therefore not produce a strong GW signal 
through the usual mechanism (see however [9]). This can however be changed in models beyond 
the SM. Extensions of the SM which lead to a strong first order electroweak phase transition are 
particularly attractive since they can provide one missing ingredient for generating the observed 
baryon asymmetry of the universe. GW signals from such models have for example been studied 
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Figure 1: Phase diagram of QCD at zero chemical potential (schematic). The dashed region 
represents our current lack of knowledge about the order of the PT in the limit of two massless 
flavours. 


in [10-16]. It is more difficult to modify the QCD phase transition, although a large neutrino 
chemical potential could be sufficient to provide a strong first order PT [17]. The resulting signal 
was studied in [18]. 

The aim of this work is to point out that gravitational waves could also be produced by a 
strong PT in a dark or hidden sector. The particular scenario we have in mind is a dark sector 
with a new SU(A^rf) gauge interaction which confines at some scale A^- Such models have recently 
received renewed interest either as models of dark matter [19-44] or as part of the low energy 
sector of so called Twin Higgs models [45-51]. Different from generic hidden sectors [52], these 
models provide a preferred mass range and some restrictions on the particle content, such that 
the frequency range of the potential GW signal can be predicted. 

Given that the SM QCD transition is not first order, we will review the known results on the 
order of the PT in strongly coupled gauge theories in the next section, followed by a discussion of 
models that fall into this category. In Sec. 3 we calculate the GW spectra that can be produced 
in these models, and compare them to the sensitivity of current and planned GW detection 
experiments in Sec. 4. We discuss the complementarity of GW experiments with other searches 
for dark sectors in Sec. 5, before presenting our conclusions. 

2 Models with First Order Phase Transition 

Near the QCD confinement scale AqcD) tbe dynamics of QCD is governed by three flavours, 
two of which are almost massless, while the strange quark mass is of order Aqcd- Lattice 
studies [5,6,53] have shown that for these values of the quark masses, the QCD PT is a weak 
cross-over. 

However this is not a generic result for QCD and similar theories, but more a consequence 
of the precise values of ~ nid and rus in the SM. The QCD phase diagram for arbitrary 
rriu^d and nis can be summarised in the so called Columbia plot, which is reproduced in Fig. 1, 
based on [54]. The pure Yang-Mills limit mu,d,^s —)• oo is known to have a strong hrst order 
PT [55] from the restoration of a global center symmetry at low temperatures. The opposite 
—)• 0 limit, i.e. the theory with three exactly massless quarks, also features a strong 
hrst order transition, related to the breakdown of the SU{3) x 517(3) chiral symmetry [56]. 
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Here we are interested in SU{N(i) theories with > 3 and Uf massless^ Dirac fermions in 
the fundamental representation, and with a confinement scale To guarantee the existence of 
a confining phase we further impose nj < 4Ai^, in order to stay outside of the conformal window. 

For pure Yang-Mills theories, n/ = 0, the confinement phase transition is related to the 
restoration of a global Z]\r^ C SU{N(i) center symmetry, which is broken in the high temperature 
phase. Lattice simulations have shown that this PT is strong first order for > 3 [57]. For 
the Nd = 3 case there is also the classic argument of [55]: In the case of a second order PT, the 
critical behaviour of the theory should be described by a Z 3 symmetric effective theory at an 
infrared fixed point. Since no fixed points are known for theories with Z 3 symmetry, the PT 
must be first order. 

The case of non-zero Uf for = 3 and —>■ 00 is discussed in [56]. For Uf = 1 there is no 
symmetry breaking, and therefore no phase transition. For nj > 3 the PT is first order^ for any 
value of Nd > 3. In the large limit, the PT is also first order for n/ = 2. The only case that 
is currently not understood is nj = 2 for finite where no analytic arguments can be applied 
and lattice results are difficult to obtain. 

To summarise, SU{Nfi > 3) theories with nj massless flavours have a first order PT if either 
n/ = 0or3<n/< 

The above discussion hopefully convinced the reader that strong PTs are generic in extensions 
of the SM which feature new confining gauge symmetries. In the following we will discuss a 
few examples which are well motivated either from a dark matter perspective or by naturalness 
arguments. These physically motivated scenarios furthermore provide a preferred range for the 
confinement scale which allows us to estimate the temperature of the PT, T* ~ A^, and 
therefore make a prediction for the GW spectrum. 

Composite Dark Matter I (CDMl) In this class of models, the dark matter candidate 
is the lightest baryonic bound state of a SU{N(i) dark sector with Uf dark quarks which are 
neutral with respect to all SM interactions. This allows the DM to be light, since it will only 
communicate with the SM through heavy mediators. The most natural realisation of these 
models is in the context of asymmetric dark matter, where the DM number density is related to 
the baryon asymmetry (see e.g. [58-60] for reviews). The measured DM density then implies a 
mass range of 0(5) GeV for the DM particle, and therefore motivates T* ~ A^ ~ (1 — 10) GeV. 

As a concrete benchmark we will consider the dark QGD model [33,43], which, as the name 
suggests, consists of a SU{Nd = 3) theory with n/ > 3 dark quarks, such that a strong PT is 
guaranteed. 

Composite Dark Matter II (CDM2) Similar to the previous case, the DM is a baryonic 
bound state of a new SU (Y^) theory. The important difference now is that the dark quarks carry 
electroweak quantum numbers. Models of these type were considered for example in [24,25], where 
mass spectra and form factors are calculated from first principles on the lattice. Constraints on the 
invisible width of the Z-boson immediately suggest a confinement scale at or above the 100 GeV 
scale, while unitarity arguments on DM freeze-out place a limit of A^^ < 0(100) TeV. For fermionic 
DM (i.e. for odd Y^), direct detection furthermore imposes a bound of Mdm ~ A^ > 10 TeV [24]. 

^More precisely, n/ is the number of fermions with m <C A^. 

^The behaviour can be nnderstood by studying the effective theory near the confinement scale, which consists of 
rij massless Goldstone bosons from the breaking of the SUinj) x SU{nf) x U{1)a chiral symmetry. For Nd —>■ 00 , 
the anomaly is negligible, and the PT is first order due to the absence of infrared fixed points in the effective 
theory, provided that nj > \/3. For finite Nd and Uf = S the anomaly contribution to the effective lagrangian (i.e. 
the rj' mass) is cubic in the Goldstone fields, and therefore alone snfficient to render the PT first order. For n/ > 3 
the anomaly does not affect the PT behaviour, so the large Nd result applies and the PT is first order. What 
remains unresolved is the rif =2 case for finite Nd- For more details the reader is referred to [56]. 
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For scalar DM (even N^) lower DM masses are allowed [25]. The benchmark model considered 
in [25] has = 4 with n/ = 4 flavours and should therefore feature a strong first order PT.^ 

Twin Higgs (TH) These models attempt to solve the hierarchy problem without introducing 
coloured partners for the top quark [45,46], and have recently received renewed attention [48-50]. 
This is achieved by adding a twin sector to the SM with an approximate Z 2 parity symmetry, 
with the minimal requirement that there should be fermionic partners for the top and bottom 
quarks which are charged under a new SU{Nd = 3) interaction. The approximate Z 2 symmetry 
constrains to lie in the (1 — 10) GeV range. Furthermore the twin top and bottom partner 
masses are necessarily much larger than A^, such that the theory at the PT scale is SU{3) with 
Uf = 0, which again predicts a strong first order PT. 

Strongly Interacting Massive Particle (SIMP) The key idea behind SIMP DM is that 
the relic density is determined through freeze-out of 3 —)• 2 instead of 2 —)• 2 annihilations [61]. 
The correct relic abundance is then obtained for DM masses in the 100 MeV range. One simple 
realisation of this mechanism is a SU{Nii) dark sector with dark pions as DM candidates. The 
3 —)• 2 processes can then be induced by the Wess-Zumino-Witten term [62,63]. The condition 
for the existence of the WZW term is nj >3, which again allows for a strong first order PT, this 
time in the 100 MeV — GeV range. 

All these models should give rise to a GW signal from the strong PT, and the resulting 
spectra and detection possibilities will be discussed in the following sections. Before moving 
on, we should mention a few other scenarios which also contain new strongly coupled sectors. 
The first group consists of Hidden Valley [52] and vector-like confinement models [64] , which 
propose the existence of new confining sectors, communicating with the SM either through heavy 
mediators or directly via SM gauge interactions. If the particle content is such that the models 
fall into one of the strong PT regions, then also these models will give rise to a GW signal. 
Hidden sectors also play a role in conformally invariant extensions of the SM [65,66], and the 
presence of a strong PT in these models was demonstrated using NJL methods in [67]. 

Furthermore there are so called mirror-world models [68], where the dark sector consists of 
an exact copy of the SM. Here one does not expect a GW signal for two reasons. First, if the 
quark content is an exact copy of QCD, there will not be a first order PT. If for some reason 
the quark masses are modified, such that a strong PT happens, cosmological constraints on the 
mirror photon require the mirror sector to have a lower temperature. Since the energy density 
is proportional to the fourth power of the temperature, and the intensity of a GW signal is 
proportional to the energy density, this will strongly suppress any GW signal from a mirror 
world. 

3 Gravitational Wave Spectra 

Gravitational waves produced at a time (or equivalently at a temperature T*) will propagate 
undisturbed in the expanding universe, therefore their frequency / and their fraction of the 
critical energy density Hqw will decrease as a~^ and a~^, respectively, where a{t) is the scale 
factor [10,69]. Denoting by a* and oq the scale factors at time of production and today, entropy 
conservation {sa^ = const) implies 


( 1 ) 

\9*,s J 

®In [25] the quark masses were taken close to the confinement scale to allow a faster simulation. In this case the 
general arguments for the order of the PT do not apply, and a lattice simulation might be necessary to determine 
the order of the PT. However the model is also viable phenomenologically with lighter quark masses. 
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Here Tq = 2.725 K = 2.348 x 10 GeV is the temperature of the CMB, igo,s = 3.91) is the 
effective number of relativistic degrees of freedom contributing to the total entropy at the time of 
production (today), and the entropy density at temperature T is given by s(T) = gs(T)T^, 
where 


2=bosons ^ ^ 2=fermions ^ ^ 


( 2 ) 


and gi counts the internal degrees of freedom of the i-th particle. It follows that the frequency 
today can be expressed as 


CL* f* 7 

/ = = 1.59 X 10"^ Hz X 

ao ri* 
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where we have used the Hubble rate at time of production, = \J j^, and assumed that 
all species are in thermal equilibrium at T = T*, i.e. 5 * = g-t^g- For the fraction of energy density 
in gravitational waves today we similarly obtain 

^GW = ^*GW = 1-77 X f —^ H*gw > ( 4 ) 

Pent \aoJ \g*J 

where we used that Pcnt/p*cnt = Hq/H^ and 77o = 2.13 x /i x 10“^^ GeV. It is noteworthy that 
the intensity of the GW signal is independent of the production temperature T* (except for the 
implicit dependence of g^ on T*). 


The most sensitive frequency regions of pulsar timing arrays and satellite based experiments 
are in the nano-Hz and milli-Hz range, respectively. To get an idea about the detectability of 
GWs from a strong dark PT we will therefore need to understand the spectrum of the produced 
GWs. For this, we will closely follow the discussion of [18]. 

Gravitational Waves are sourced by tensor fluctuations of the energy momentum tensor of 
the primordial plasma. During first order phase transitions both bubble collisions [70-72] and 
magnetohydrodynamical (MHD) turbulence [73-79] provide sources of gravitational waves. As 
functions of the conformal wave number k = 27ra/, the GW spectra produced by either source 
can be approximated by [18] 


h? 


dlog k 
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1 + (fc//3) 
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{k/l3f 


(1 -I- Ak/l-L^) (1 -I- {v/Tt‘^){k/ 


( 5 ) 

( 6 ) 


Eqn. (5) is based on [80,81] while Eqn. ( 6 ) is a fit to the numerical results of [82]. Here Td* is 
the conformal Hubble parameter T-L = Ha at T = T*, and is the radiation energy density 
today. The quantities that determine the GW spectrum are the bubble nucleation rate /3 (the 
duration of the PT is the bubble velocity v and the relative energy density in the source, 

^S* = Ps*/p*,crit = The temperature of the PT enters through the dependence of H* on 

H. 

The duration of the PT is usually taken as (1 — 100)% of a Hubble time, and therefore 
/3 = (1 — 100)77 [2]. To understand the relation with the physical frequency, remember that the 
conformal frequency is related to the conformal wave number via af = k/{2Ti). Eurthermore 
using H = Ha we see that f*/H^ = af^/H* = (k/H*)/{27r), which together with Eqn. (3) allows 
us to translate the GW spectra into physical frequencies. 

In a given theory, the dynamics of the phase transition, and therefore the parameters (3, v 
and 0 , 3 *, are in principle calculable. Eor the strongly coupled models considered here they are 
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Figure 2: Left: Peak frequencies of the GW spectra (in mHz) from bubble collisions (blue,solid) 
and MHD turbulence (red, dashed) in the T* — /? plane, for v = 1.0. Right: GW spectrum from 
bubble collisions (blue, solid) and turbulence (red, dashed) as well as the combined spectrum 
(black, thick), as a function of conformal wave number k, for v = 1.0 and /3 = lOT-L^. 


however not known, and can only be estimated using lattice simulations. We will therefore take 
/?, V and ils* as additional input parameters, with values motivated by results of analyses in 
weakly coupled models. 

Following [18], we will use Qs* = 0.1 and j3 = lOL^*, but v = 1.0 instead of 0.7. We are now in 
a position to study the location of the peaks of the GW signals from bubble collisions and MHD 
turbulence. The bubble collision signal is triangular shaped with a maximum at A://3 = -^3 ~ 1.3, 
while the MHD turbulence peaks at somewhat larger wave numbers A://3 ~ vr^/u. To obtain 
physical frequencies, we use Eqn. (3) and = {l3/{2'ir'H)){k/f3). Then the peak locations are 


/S = 3 33 X 10-» Hz X 




AMHD) _ AB) 

J peak ~ peak * 


( 7 ) 


In Fig. 2 we show the location of the frequency peaks as function of the PT temperature T* 
and /3. As expected from Eqn. (3), the peak frequencies increase linearly with the transition 
temperature T* and with /l/H*. 

The source term D 5 * can be different for bubble collision and turbulence. Here we will assume 
that equal amounts of energy act as source for and . In this case the turbulence 

signal dominates over the one from bubble collisions over most of the relevant frequency range, 
see Eig. 2. The intensity of both signals decreases as (/3/?7*)“^, therefore smaller values of j3 
are preferable. From Eqn. ( 6 ) it might appear that the turbulence signal only decreases as 
however the k/T-L^ term in the denominator gives another power of for k^l. 


Recent simulations of first order PTs suggest that sound waves generated by the expansion 
of bubbles could be the dominant source of GWs from these transitions [83-85]. Sound waves 
continue propagating through the early universe after the PT is finished, and decay on a timescale 
"R*. Gompared to the above discussed spectra, they will therefore not be suppressed as much by 
the velocity of the transition /3, and the signal could be increased by a factor (/3/77*) compared 
to the bubble collision signal, but with a spectrum decaying as k~^. This could potentially boost 
the signal, in particular for cases where the PT is fast, i.e. /3/'H 1. 
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Figure 3: GW spectra n{f)h^ for T* = 0.1 GeV (SIMP), T* = 3 GeV (GDMl, TH models), 
r* = 300 GeV and T* = 10 TeV (GDM2 models). The upper (lower) edges of the contours 
correspond to (3 = % (/? = lO'H), and furthermore v = 1 and fig* = 0.1 for all curves. The red 
band T* = 0.1 GeV indicates where a signal of the QGD PT would lie if it was strong. The 
projected reach of several planned GW detection experiments is shown (dashed). 

4 Detectability 

In the previous section, we have seen that the peak frequencies of GW signals from GeV-TeV 
scale PTs are of order (10“® — 10“^) Hz. Furthermore it is important to note that a broad 
spectral region around the peak is populated by GWs, from (10“^^ — 1) Hz. 

GWs with frequencies down to 10“® Hz can be probed by satellite based experiments like 
eLISA [86], however the sensitivity quickly degrades below 10“^ Hz. On the other end of the 
spectrum, pulsar timing arrays (PTA) can probe frequencies in the (10“® — 10“^) Hz range. In 
Fig. 3 we overlay the expected GW signal for different model parameters with the expected 
sensitivities of current and planned GW detection experiments (based on [87]). 

Glearly the most promising signals are those from models with a PT temperature in the TeV 
range, where the peak region of the GW spectrum falls right into the most sensitive frequency 
range of satellite experiments. Here the signal should even be detectable for choices of the 
parameters that are less optimistic than those used for Fig. 3. Models of the GDM2 type 
naturally fall into this region, but also the GDMl models can be viable with a conhnement scale 
in that region. 

The fact that TeV dark sectors predict an observable GW signal is not surprising, since here 
we are in the energy range of the electroweak PT or beyond, and the observability in particular 
of TeV scale strong PTs has been noted before [89-91]. The novelty here is that the dynamics 
leading to this strong PT does not have to be connected to the electroweak sector of the SM, 
and is therefore not in tension with the non-observation of new physics at the LHG. 

For models with T* ~ (1 — 10) GeV the situation is a bit more difficult, since the signal peak 
ends up in a frequency region where neither PTA nor (e)LISA are sensitive. Looking at the 
r* = 3 GeV curve in Fig. 3 more closely, we see that in the best case scenario, for /3 = 77, both 
PTAs and LISA would be able to detect parts of the GW spectrum. For larger f3 the signal 
quickly drops out of the PTA sensitivity region, however LISA remains sensitive. This is due 
to the increase of the observed frequency with which partially compensates the overall 

(77//3)^ drop of the signal in the LISA sensitivity region. Therefore there is a chance to detect a 
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GW signal from the CDMl and TH models, even if the predicted frequency range is not optimal. 

Recently the EPTA experiment has reported the first limit on a stochastic GW background [88] , 
which probes values of of order 10“^ in the (10“® — 10“^) Hz range. Their limit can not 

directly be displayed in Fig. 3 since it depends on the assumed spectrum, but it is nevertheless 
interesting to note that GWs from a 100 MeV scale dark sector could already be detectable. 

Another way of interpreting these results is that GW searches have the potential to observe 
signals of completely unknown physics in the early universe, while the absence of a signal in a 
certain frequency range would place (very weak) bounds on the dynamics of new physics at the 
corresponding energy scale. 

5 Complementarity 

Signals of the new physics models introduced in Sec. 2 are being searched for at collider 
experiments and in direct and indirect dark matter searches. The possibility to observe such 
a signal in those experiments always relies on sufficiently strong non-gravitational interactions 
of the dark sector with the SM. Instead the GW signal is unique in the sense that it probes 
the gravitational effects of a model at very early times, and, while the detectability of a signal 
eventually still depends on model parameters, those are relatively independent of whether the 
model is detectable at colliders or in dark matter searches. 

In GDMI type models, or more general, hidden valley models with a new strong interaction, 
the dark sector is neutral under all SM interactions, and communicates only through a heavy 
mediator. If that mediator is in the TeV range, both direct detection experiments and the LHG 
could discover these models, but for larger masses it becomes increasingly difficult. From the 
DM perspective on the other hand, mediator masses of 10s of TeV are acceptable, such that part 
of the parameter space will remain unexplored in the near future. The GW signal instead is 
completely independent of the mediator mass, such that it could be detected even if the model 
was not discovered before, and thanks to the straightforward connection of mass scales and GW 
frequencies, could even motivate future experimental efforts. 

The situation is similar for TH models, which mainly communicate with the SM through 
the Higgs portal. Here the parameter that controls detectability is not a mass scale but the 
smallness of the mixing angle of the Higgs with its twin partner. First hints for such a model 
could come from deviations in Higgs couplings and from exotic Higgs decays, but even then it 
would be difficult to uncover the whole structure. A GW signal in the right energy range could 
provide much information about the dynamics of these models near their confinement scale. 

Finally the GDM2 type models have many detectable features at hadron colliders. However 
the upper limit on their mass scale is of order 100 TeV, which comes from the requirement 
that dark matter is not overproduced. Such high scales are not in reach of the current or next 
generation of collider experiments, but seem very accessible by GW searches, since the signal 
moves into the most sensitive frequency region of LISA and other satellite based experiments. 

Overall we see that GW experiments provide a unique window to explore the dynamics of 
these models in the early universe, even if they are not discovered in the near future. Of course 
it would be even more interesting if such a dark sector was discovered at the LHG. The GW 
signal would then provide an important and unique probe of the physics of those models in the 
early universe. 

Moving beyond the concrete models discussed above, perturbative unitarity constrains the 
mass of thermal DM to be below 110 TeV [92,93]. For composite non-perturbative DM this 
limit does not apply directly^, instead a lower bound on the radius of the extended object can be 
obtained, R < (100 TeV)“^. It is reasonable to expect the radius R to be of order of the inverse 
mass, which again implies an upper bound on the DM mass of order 100 TeV. GW signals could 
therefore be a unique probe of the thermal DM paradigm. 

^We thank T. Cohen for discussion of this point. 
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6 Conclusions 


Models beyond the standard model with a confining dark sector can lead to unexpected phe¬ 
nomenological signatures. Here we have explored the possibility to detect gravitational waves 
due to a first order phase transition at the confinement scale A^- The main messages from this 
paper are as follows: 

• Different from QCD, dark sectors with QCD-like interactions can undergo strong first order 
phase transitions, with only mild constraints on the particle content of the theories. 

• Several classes of new physics models that are currently being explored fulfil the criteria 
for first order PTs. The physics problems these models are trying to address, either dark 
matter or naturalness, constrain the confinement scales and therefore the temperature 
range of the phase transition. 

• The GW signals originating from these dark phase transitions are in the detectable frequency 
range of future GW experiments, either at (E)LISA and BBO for high scale models, or in 
PTA experiments for the lower end of the spectrum. 

Depending on other aspects of the model, GW signals will either provide complementary 
information about the models in question, or might even be the the best option to find evidence 
for these models of new physics. Different from the electroweak PT, a null result at LHC will 
not strongly disfavour a strong PT in a dark sector, although of course a confirmation of their 
existence would be more exciting. 

A shortcoming of the present work is a lack of precise quantitative predictions. The bubble 
velocity v as well as the time scale of the phase transition /3“^ and the energy fraction are 
currently unknown, and are set to optimistic (but not unrealistic) values. Two approaches seem 
possible to improve upon this situation: On one side, lattice simulations could be used to measure 
quantities like the latent heat and the surface tension, which are related to the above parameters 
and can be used to obtain a more quantitative prediction for the GW spectra. Alternatively, one 
could attempt to construct a holographic dual for some of these theories, and analyse the PT in 
that setup. 
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